The CH3SH molecule deposited on Cu(lll) and deprotonation 
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Abstract 

We demonstrate for the first time that when a methanethiol adsorbed on the regular Cu(lll) 
surface, the dissociative structure is thermodynamically more stable than the intact one. The 
computational results show that at low temperature the methanethiol adsorbate prefers the atop 
site of the regular Cu(lll) surface. As the temperature is increased, the S-H bond is broken 
and the methylthiolate favors the hollow sites. On the defected Cu(lll) surface, the dissociative 
configuration is still thermodynamically more stable than the nondissociative one. The calculation 
indicates that the hydrogen initially attached to the sulfur would like to form a bond with the 
copper surface rather than desorb from it. Even though both copper and gold are the noble metal, 
the stability of the methanethiol adsorption on the Cu(lll) substrate is almost the reverse of that 
on the Au(lll). 

PACS numbers: 61.46.-w, 36.40. Cg, 68.43.Bc, 73.20.Hb 
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I. INTRODUCTION 



In recent years, a lot of work has been conducted on the interaction of alkanethiols with 
gold surface because of their ability to form self-assembled monolayers (SAMs); see Ref. 
l| for a review. These SAM films have wide application on wetting phenomena, tribology, 
chemical and biological sensing, optics and nanotechnology. The consistent model for the 
methanethiol adsorption process on the Au(lll) substrate was proposed recently RQ. It 
was revealed that the methanethiols stay intact on the regular Au(lll) surface, but the S-H 
bond ruptures on the defected Au(lll) surface. On other noble metals such as copper, the 
adsorption pattern of the alkanethiols (or alkanethilates) on the Cu(lll) is distinct from that 
on the Au(lll) surface {4]- 15|. Then one may ask how the S-H bond in the methanethiol 



on the regular (or defected) Cu(lll) surface breaks. 

The normal incidence X-ray standing wave measurements show that at the low tem- 
perature, the intact methanethiols prefer a local atop site on the Cu(lll) surface, but at 
140K methanethiols turn into methanethiolates depositing on unreconstructed areas of the 
Cu(lll) surface, occupying two nonequivalent hollow sites with similar probability IJ]. At 
room temperature, when methanethiols or dimethyl disulfides adsorbed on the Cu(lll), 
the formation of the pseudosquare reconstruction was observed via scanning tunneling mi- 
co^copy BQ. A grazing mcMence surface X-.ay study iud.cate. that intact alkanethiol. 
favor hollow sites on Cu(lll) [7|. However, there has been no experimental proof for hydro- 
gen desorption from such processes. As we know, alkanethiols only stay intact on the regular 
Au(lll) surface, but the S-H bond ruptures on the defected Au(lll) surface . So far, 
no theoretical work on the scission of the methanethiol on the Cu(lll) surface has been 
done, thus it is imperative to study the mechanism of the S-H bond rupture theoretically. 

In this paper, we study the adsorption of methanethiol on the Cu(lll) surface by the 
density functional theory 16|]. First, we will present adsorption energies and geometries 
for the nondissociative and dissociative configuration on the regular Cu(lll) surface. We 
demonstrate, for the first time, the dissociative structure is the most stable one for the 
regular Cu(lll) surface ^hich is opposite to the case when the methanethiols adsorbed on 
the regular Au(lll) j^, 3]. The S atom of the intact methanethiol stays on the top of the 
copper atom of the surface. The energy barrier is about 0.2 eV. This implies that at low 
temperature, the intact methanethiol prefers a local atop site on the regular Cu(lll) surface 
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41. As the temperature increased, the energy barrier can be overcome and the hydrogen is 

nriQ 

detached from the S atom, which is consistent with the experimental observation |J, 16|, |7|. 
To discuss the stabihty of configurations for the defected Cu(lll) surface, we assume a 
vacancy in the top layer of the Cu(lll) surface. We find that on the defected surface, the 
dissociative structure is still more stable than nondissociative one. The hydrogen initially 
attached to the sulfur forms a bond with the Cu(lll) substrate instead of desorbing from 
it. 



II. COMPUTATIONAL METHOD 



The calculations were performed by the VASP code [l^ in two approaches: 1) using the 
projector augmented wave potential |l7|, [l8| and the Perdew-Wang 91 (PW91) generalized 



gradient approximation to exchange and correlation (XC) [19']; 2) the ultrasoft pseudopo- 
tential [l^ with XC based on the Perdew-Burke-Ernzerhof (PBE) formalism [20] . The wave 
functions are expanded in a plane wave basis with an energy cutoff of 400 eV. The Brillouin 
zone integration is performed by use of the Monkhorst-Pack scheme 21j. We utilized a 



3 X 3 X 1 k point mesh for the geometry optimization. The adopted supercell corresponds 
to 12 Cu atoms per layer. The Cu atoms in the top three atomic layers are allowed to relax, 
while those in the bottom layer are fixed to simulate bulk-like termination 22i]. The vac- 
uum region comprises seven atomic layers, which exceeds substantially the extension of the 
methanethiol molecule. To examine the accuracy of our approach we increased the energy 
cutoff to 500 eV and the number of k points to 8 x 8 x 1. In both cases, the difference 
amounted to less than 1.8%. We further calculated the copper lattice constant, and found 



it to agree with the experimental value 



23j within 0.6%. 



III. RESULTS AND DISCUSSION 



A. Methanethiol molecule deposited on the regular Cu(lll) surface 

We begin with the geometries and adsorption energies of the optimized structures for 
the methanethiol on the Cu(lll) surface at 0.25 ML, as displayed in Table [B The symbol 
fcc-bri (or hcp-bri) in Table [T] represents that the S atom is placed in a fee (or hep) hollow 
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TABLE I: The geometries and adsorption energies for the optimized nondissociative and disso- 
ciative configurations on the regular and defected Cu(lll) at 0.25 ML. The entries S site, 9, tilt, 
ds-Cu (^) and Eads (eV) refer to the S headgroup atom site, the angle between the S-C bond 
direction and the normal to the Cu(lll), the region of th S-C bond tilted, the shortest S-Cu bond 
length and the adsorption energy, respectively. 
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FIG. 1: (color online), (a) Intact structure on the regular Cu(lll). (b) Dissociative structure in 
the fee hollow site of the regular Cu(lll) (fcc-bri). (c) Dissociative structure in the hep hollow site 
(hcp-bri). (d) Dissociative structure on the Cu(lll) with a vacancy. 

region, but leaned toward the bridge. The stable configurations on the Cu(lll) surface are 
shown in Fig. [H 

We calculated nineteen nondissociative configurations on the regular Cu(lll) surface, 
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which are characterized by the S atom adsorption site (atop, bridge, fee and hep hollow), the 
angle of the S-C bond, and the direction of the H-S bond. The energy difference among the 
configurations is around 0.1 eV, i.e., the adsorption energies span 0.30 eV - 0.61 eV. However, 
the adsorption energies weakly depend on the tilting of the alkyl chain, which reveals the 
ionic nature of the bonding. The first row of Table [I] displays that the adsorption energy 
for the most stable nondissociative structure is 0.61 eV (Fig. [1^), and the S adsorption 
site is on the top of the Cu atom. Thus at low temperature, the intact methanethiol 
stays on the top site of the regular Cu(lll) surface, in keeping with a recent experimental 
result ^. The second and third rows of Table [T] list two most stable dissociative structures 
on the regular Cu(lll) substrate, and their adsorption energies are 1.14 eV and 1.11 eV. 
The S atom for both configurations are in hollow site, i.e^ fcc-bri (Fig. [TJo) and hcp-bri 
(Fig. [It), which is in favor of the explanation in Ref. Il2| (when dynamical effects are 
considered, only the hollow site occupancy is observed [12;]). The first and second rows 
in Table [T] reveal that the adsorption energy of the dissociative structure is higher than 
that of the nondissociative one by 0.53 eV. To further examine this result, we apply the 
ultrasoft pseudopotential in conjunction with the PBE XC j2o| to recalculate the above 
structures, and find an adsorption energy difference of 0.51 eV between the cases of cleavage 
and noncleavage. Thus we conclude that the dissociative configuration is more stable than 
the nondissociative one, which is consistent with the experimental observation that at low 
temperature, the intact methanethiol stays on the regular Cu(lll) surface, and when the 
temperature increased, methanethiolates deposited on unreconstructed areas^of the Cu(lll) 
surface, occupying two nonequivalent hollow sites with similar probability 



B. Projected density of states 

To see why the dissociative structure is more stable than the nondissociative one, we 
calculated the partial density of states (PDOS) for the H atom initially attached to the 
sulfur atom (Fig. [2^). The PDOS projected on the S atom adsorbed on the regular Cu(lll) 
surface is illustrated in Fig. [2]d. Fig. [2^ displays that the primary peaks of the hydrogen 
PDOS for the dissociative structure are at positions of lower energy than those for the 
nondissociative one. The electronic configuration related to the former case is thus more 
stable than that corresponding the latter, and the hydrogen atom forms a stronger bond 
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FIG. 2: (color online), (a) The PDOS projected on the hydrogen atom initially attached to the S 
atom adsorbed on the regular Cu(lll) surface, (b) The PDOS projected on the S atom adsorbed 
on the regular Cu(lll) surface. 

with the Cu(lll) surface than with the sulfur atom. 



C. Interacting bond 

To elucidate the interacting bond between the methanethiol (or methylthiolate) and the 
Cu(lll) substrate, we evaluate the spatially resolved charge-density difference: Ap(r) = 
Pads/sub{r) - Psub{r) - Pads{r), where Pads/sub, Psub and pads are the electron charge densities 
of the relaxed adsorbate-substrate system, of the clean relaxed surface, and of the isolated 
but adsorption like deformed adsorbate (without substrate), respectively. The isodensity 
surfaces of the charge-density difference for the intact and dissociative structure are depicted 
in Fig. Fig. [3)3 and Fig. [2t (to discuss the change of the S-Cu bond before and after 
dissociation, in Fig. [3l we display only the surrounding part of the S-Cu bond). 

The calculated charge transfers from the Cu(lll) substrate to the methanethiol for the 
nondissociative adsorption and to the methylthiolate for the dissociative structure are 0.3e 
and 0.6e, respectively. Fig. [3^ shows that for the intact adsorption, the electrostatic inter- 
action responsible for the bonding comes from the mono pole term (the 0.3e charge transfer) 



and the dipole moments in the adsorbate and substrate 2J] . Around the S atom, there is a 



"ring" of accumulation of electron charge. The electrostatic interaction is dominated by the 



attractive ionic term modified by a repulsive dipolar term 2J]. The sulfur stays on the top 
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FIG. 3: (color online). The isosurfaces of the charge-density difference for (a) nondissociative ad- 
sorption with blue (accumulation of electrons) /red (depletion of electrons) isosurface value: ±0.015 
A , (b) dissociative structure (front view) and (c) dissociative structure (back view) with isosur- 
face value: ±0.025 e/A . Only three related copper atoms of the Cu(lll) surface are illustrated. 



of the copper atom, i.e., the sulfur only forms a bond with one Cu atom. The interaction 
between the methylthiolate and Cu(lll) is the ionic bonding 10|, ll5|. The sulfur in the 
methylthiolate forms bonds with two copper atoms of the Cu(lll) surface (see Fig. [3]d and 
Fig. [St), and its charge transfer (0.6e) is two times of that for the nondissociative case 
(0.3e). It indicates that the methylthiolate takes in 0.3e electron charge from each copper 
atom channel (0.3e x 2 = 0.6e), which matches with the fact that the adsorption energy for 
the dissociative structure (1.14eV) is about two times of that for intact adsorption (0.61eV). 
The charge densities in the [1-21] and [10-1] plane for intact and dissociative adsorption are 
shown in Fig. |U which is consistent with the picture obtained from the isosurfaces of the 
charge- density difference, charge transfers, and bond structures. 



D. Methanethiol molecule deposited on the defected Cu(lll) surface 

On the regular Cu(lll) surface, we found that the dissociative structure is thermody- 
namically more stable than the nondissociative one. To see if the defects can change the 
above conclusion, we consider the Cu(lll) surface with a vacancy in the top layer of the 
supercell. The geometries and adsorption energies for the optimized nondissociative and 
dissociative configurations on the Cu(lll) surface with a vacancy are indicated in the fifth 
and sixth rows of Table [B It shows that the adsorption energy for the dissociative structure 
(Fig. [T]i) is 0.9 eV higher than that for the intact one, i.e., the dissociative configuration 
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[10-1] 



[1-21] 




FIG. 4: (color online). The contour plots of the charge densities, (a) Charge density plot in 
(1-21) plane which extends along [111] and [10-1] directions and (b) in (10-1) plane for the intact 
adsorption, (c) Charge density plot in (1-21) plane and (d) in (10-1) plane for the dissociative 
adsorption. 

on the defected Cu(lll) surface is thermodynamically more stable than the nondissociative 
one. For the stable dissociative adsorption, the S atom forms two bonds with two adjacent 
copper atoms which surround the vacancy. 
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Transition State 
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FIG. 5: (color online). The structures for the initial intact adsorption, the final dissociated one, 
and the six images. 

E. Energy barrier 

To explain why at low temperature, the intact methanethiol on the regular surface was 
observed, and around 140K, the S-H bond in the methanethiol on the Cu(lll) was broken, we 
calculated the energy barrier for the process that lead from the nondissociative structure Fig. 
[T^ to the dissociative structure Fig. [Tb. The energy barrier was calculated by the climbing 



nudged elastic band method, [25|-[27|] where six equidistant images have been used. The 
atomic configurations for the initial nondissociative adsorption, the six intermediate images, 
and the final dissociative one are shown in Fig. [51 The structure of the transition state is 
depicted in Fig. [Sji. 

The relative energies along the S-H bond cleavage reaction path are schematically illus- 
trated in Fig. O The energy barrier is about 0.2 eV, which is lower than the corresponding 
energy barrier for the scission of the S-H bond in a single methanethiol molecule. The rea- 
son is that when the methanethiol is adsorbed on the Cu(lll), the surface plays a role as a 
catalyst, which can suppress the barrier height. As the temperature is raised, the S-H bond 
break reaction is activated over the energy barrier, which is consistent with the experimental 



measurement 



a. 
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FIG. 6: The relative energy levels along the S-H bond scission reaction paths. 
F. Hydrogen atom behavior 

We have seen that at low temperature, the intact methanethiol favors the atop site of the 
Cu(lll) surface, and when the temperature increased, the S-H bond in the methanethiol on 
the Cu(lll) surface is ruptured. The hydrogen could either stick to the Cu(lll) surface or 
desorb in molecular form. No experimental evidence of either process has been found. To 
predict if the hydrogen forms a bond with the Cu(lll) surface or desorbs in the hydrogen 
molecule, we compared the interaction energy for the dissociative structure E(^ch3S+h)/Cu(iii) 
with that for methylthiolate on the Cu(lll) and the desorbed hydrogen molecule, i.e., 
E(CH3S)/Cu{iii) + ^Eh2- We found that E(_ch3S+h)/Cu{iii) is lower than E(^ch3S)/Cu{iii) + 
\Efj^ by 0.24 eV, which implies that the hydrogen would like to form a bond with copper 
surface rather than desorb from the surface. 

To check if the deposited hydrogen modified the methylthiolate adsorption pattern, we 
calculated fifteen methylthiolate structures on the regular Cu(lll) surface. The most stable 
structure is shown in the fourth row of Table [H The second and fourth rows in the Table 
display that the S-Cu bond lengths and the angles 9 for the dissociative structure (both 
the hydrogen and the CH3S adsorbed on the Cu(lll) surface) and for the methylthiolate 
(only CH3S) on the Cu(lll) are 2.26 A, 40.2°, 2.28 A and 46.2°, respectively. The hydrogen 
adsorbed to Cu(lll) modifies the methylthiolate adsorption structure by 6°. 
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The regular Au(1 1 1 ) surface 



H attached to S 




FIG. 7: (color online). The PDOS projected on the hydrogen atom initially attached to the S atom 
adsorbed on the regular Au(lll) surface. 

G. Comparison with the methanethiol adsorption on the Au(lll) surface 



We compare the methanethiol adsorption process on the Cu(lll) substrate with that on 
the other noble metal surface - Au(lll). On the regular Cu(lll) surface, the methanethiol 
dissociative adsorption is more stable than the intact one. The corresponding energy barrier 
is about 0.2 eV. So at low temperature, the methanethiol stays intact on the atop site 
of the regular Cu(lll). When the temperature is raised, this energy barrier is overcome, 
and the S-H bond breaks. The corresponding hydrogen is adsorbed on the surface. The 
presence of the defects on the Cu(lll) does not change the situation. However, on the 
regular Au(lll) surface, the nondissociative structure is more stable than the dissociative 
one. At low temperature, the methanethiol favors the atop site of the regular Au(lll). As 
the temperature is increased, the methanethiol molecule desorbs from the regular Au(lll). 
For the defected Au(lll), this order of stabilities is reversed. To unveil the origin of the 
difference between the adsorption on the Cu and Au surface, we calculated the partial 
density of states (PDOS) for the H atom initially attached to the sulfur atom adsorbed on 
the regular Au(lll) surface [see Fig. [7j. 

Comparing Fig. [7] with Fig. [2^, we have found that for the Cu surface the primary peaks 
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of the hydrogen PDOS for the dissociative structure are at positions of lower energy than 
those for the nondissociative one. However, for the Au surface, the primary peaks of the 
hydrogen PDOS for the intact structure are at positions of lower energy than those for the 
dissociative one. Thus such a comparison explained the physical reason of the difference 
between the Cu and Au surface. 



IV. SUMMARY 



In conclusion, we have demonstrated for the first time that when a methanethiol adsorbed 
on the regular Cu(lll) surface, the dissociative structure is thermodynamically more stable 
than the intact one. The energy barrier for the process that lead from the nondissociative 
structure to the dissociative one is about 0.2 eV. Thus at low temperature, the methanethiol 
adsorbate prefers the atop site of the regular Cu(lll) surface As the temperature rises, 
the S-H bond is broken, and the methylthiolate favors the hollow sites ^. On the defected 
Cu(lll) surface, the dissociative configuration is still thermodynamically more stable than 
the nondissociative one. The calculation shows that the hydrogen initially attached to the 
sulfur would like to form a bond with copper surface rather than desorb from it. Even though 
both copper and gold are the noble metal, the stability of the methanethiol adsorption on 
the Cu(lll) substrate is almost the reverse of that on the Au(lll). 
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